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Abstract 


The  purpose  o£  this  study  was  to  evaluate  the 
performance  and  stability  robustness  of  a  lateral  autopilot 
designed  for  a  transport  aircraft  using  H2  and  Hoo  control 
theory.  The  Intent  was  to  design  a  controller  that  met 
performance  and  robustness  specifications  over  a  range  of 
flight  conditions.  First,  a  control  structure  to  be  used  In 
designing  the  autopilot  was  developed.  Once  this 
was  accomplished.  It  was  formulated  Into  the  small  gain 
problem.  Controllers  were  then  developed  using  H2  and  Hoo 
control  theory.  The  final  task  involved  evaluating  the 
controllers  developed  In  a  closed  loop  simulation. 
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ROBUSTNESS  EVALUATION  OF  H2  AND  Hoo  CONTROL  THEORY 
AS  APPLIED  TO  A  TRANSPORT  AIRCRAFT 

I .  Introduction 

1.1  Background 

In  June  of  1989,  the  results  of  a  design  challenge  were 
presented  at  the  American  Control  Conference ( 1 : 565-620 ) .  The 
focus  of  the  challenge  was  to  design  a  lateral  autopilot  for 
a  transport  aircraft.  Several  papers  were  presented 
utilizing  different  design  techniques.  These  included 
classical  control,  a  model-matching  approach,  LQG/LTR,  and  Hco 
control  theory. 

The  state  space  equations  for  fourteen  flight  conditions 
vere  provided.  In  addition,  a  list  of  design  requirements 
were  given. 

One  objective  of  the  design  challenge  vas  to  design  a 
controller  that  would  satisfy  the  design  requirements  for  all 
fourteen  flight  conditions  without  the  need  for  gain 
scheduling.  Many  of  the  autopilots  In  aircraft  and  missiles 
today  were  designed  using  single-input  single-output  (SISO) 
design  techniques  at  specific  flight  conditions.  This 
requires  a  large  number  of  gains  for  the  flight  envelope. 
Design  techniques  such  as  regulator  theory  that  handle 
multiple-input  multiple-output  (MINO)  systems  lack  the 
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robustness  to  handle  parameter  variations  and  sensor  noise. 

The  H2  control  theory  and  the  more  recent  Hoo  control 
theory  are  designed  to  be  robust  In  the  face  of  design 
uncertainties.  H2  control  theory  handles  structured 
uncertainties  such  as  parameter  variations  effectively.  The 
Hod  control  theory  Is  very  effective  against  unstructured 
uncertainties  which  Include  unmodelled  dynamics  and  sensor 
noise . 

The  above  features  of  H2  and  Hoo  control  theories  yield 
controllers  that  are  robust.  This  allows  controllers  to  be 
designed  that  reduce  or  eliminate  the  need  for  gain  scheduling 
over  a  given  range  of  flight  conditions. 

1.2  Purpose 

The  purpose  of  this  thesis  is  to  use  H2  and  Heo  control 
theory  to  design  a  lateral  autopilot  for  a  transport  aircraft 
using  the  data  from  the  design  challenge  at  the  1989  American 
Control  Conference.  The  state  space  representations  for  the 
fourteen  flight  conditions  will  be  used  in  an  attempt  to 
develop  a  controller  that  will  satisfy  the  given  design 
requirements.  The  problem  will  be  formulated  into  the  small 
gain  or  standard  problem  as  defined  in  Doyle  et.  al.  (2:631). 

A  detailed  analysis  of  the  problem  setup  and  formulation  will 
be  covered. 

The  controllers  developed  will  be  tested  to  see  if  they 
are  robust  enough  to  handle  the  parameter  variations  of  the 
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fourteen  flight  conditions.  The  H2  and  Hao  controllers  will 
be  compared . 

1.3  Overview 

In  Chapter  2,  aspects  of  multivariable  control  theory 
are  discussed  as  well  as  H2  and  Hoo  theory.  Chapter  3  covers 
the  problem  description  and  characterization  of  the  open-loop 
data.  The  formulation  of  the  standard  problem  is  contained 
in  Chapter  4.  In  Chapter  5,  the  procedure  and  analysis 
techniques  required  to  develop  the  H2  and  Hoo  controllers  are 
discussed.  In  Chapter  6,  the  results  obtained  from  running 
the  controllers  in  a  closed  loop  simulation  are  presented. 
Chapter  7  contains  a  summary  and  conclusions  followed  by 
recommendations  for  further  study. 

The  software  package  PC-MATLAB  was  used  for  this  thesis. 
This  included  the  basic  program  along  with  the  Control  System 


Toolbox . 


II.  HI  ail4  Uoo  Control  Theory 


2.1  Singular  Values  and.  ShagJjpq 

Before  H2  and  Ho  control  theory  are  covered,  it  is 
important  to  discuss  a  few  fundamentals  of  multivariable 
control  design.  One  of  these  is  singular  values. 

Singular  values  of  a  matrix  A  are  the  non-negative 

* 

square  roots  of  the  eigenvalues  of  A  A  as  given  by 

g.(A)  =  tX.  (A*A)  1^''^  (2.1) 

V  V 

iQxn  ^ 

where  AeC  and  A  is  the  complex  conjugate  transpose  of  the 
matrix  A.  The  largest  singular  value,  o^ouc,  is  denoted  ct. 

The  smallest  singular  value,  o^in,  is  denoted  a. 

For  a  transfer  function  G(j(«>),  ^[G]  and  alG]  can  be 
computed.  These  two  quantities  are  frequency  dependent  and 
give  a  measure  of  the  size  of  the  transfer  function 
matrix.  In  SISO  systems,  the  Bode  magnitude  plot 
shows  the  magnitude  of  a  transfer  function  as  a  function  of 
frequency.  The  singular  value  plot  is  the  MIMO  analog  of  the 
Bode  magnitude  plot  (4:4-11). 

The  idea  of  loop  shaping  is  also  Important.  Figure  2.1 
shows  a  typical  MIMO  system.  The  concepts  of  a  good  loop 
shape  carry  over  to  the  MIMO  system.  In  SISO  systems,  to 
have  good  command  following  and  disturbance  rejection,  the 
loop  gain  must  be  large  at  low  frequency.  A  small  loop  gain 
is  required  for  attenuation  of  modelling  errors  at  high 
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K(8) 


G(s) 


Y(s) 


Fig.  2.1  MIMO  System  Block  Diagram 
frequency.  To  avoid  sensor  noise,  which  Is  a  high  frequency 
phenomenon,  the  loop  gain  must  be  small.  Therefore,  a  good 
loop  shape  Is  as  shown  In  Figure  2.2.  A  detailed  discussion 
of  loop  shaping  Is  covered  by  Rldgely  and  Banda  (4:4-1:4-10) 


|g(  ja>)k(  jw)  I 


Fig.  2.2  Good  Loop  Shape  (SISO) 
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As  mentioned,  the  singular  value  plot  is  the  equivalent 
of  the  Bode  magnitude  plot.  The  ef[GK]  and  ^[GK]  give 
information  concerning  the  loop  gain  and  bandwidth  of  the 
system.  The  loop  shape  of  Figure  2.2  for  the  SI SO  system 
equates  to  Figure  2.3  for  the  MIMO  system.  At  low  frequency 
o'lGK]  must  be  large  to  ensure  performance  and  9[GK]  must  be 


Fig.  2.3  Good  Loop  Shape  (MIMO) 
small  at  high  frequency  to  give  desirable  robustness. 

These  ideas  will  be  used  in  designing  the  controller  for 
the  lateral  autopilot. 

2.2  112.  And.  ilo»  Control  XbfifiLX 

The  standard  problem,  as  discussed  in  Francis  (3:Chap  3) 
is  shown  by  Figure  2.4.  All  signals  can  be  vector 
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Fig.  2.4  Standard  Problem  Block  Diagram 
quantities.  The  exogenous  Inputs,  d,  to  the  plant,  G,  are 
tracking  signals  or  disturbances.  The  regulated  output  of 
the  system  Is  denoted  e.  The  measured  output  fed  back 
through  the  controller,  K,  Is  y.  The  control  inputs  to  the 
system  are  denoted  u.  The  plant  can  be  partitioned  as 


Gil  G12 
_  G21  G22 


vhlch  yields  the  equations 

e  «  Gild  *  Gl2u  (2.3) 

y  -  G21d  -t-  G22u  (2.4) 

u  »  Ky  (2.5) 

The  equivalent  state  space  form  is 

«  Ax  <!■  Bid  *  B2u  (2.6) 

e  »  Clx  *  Dlld  4  D12u  (2.7) 

y  «  C2x  *  D21d  D22u  (2.8) 


Given  the  above,  the  standard  problem  is  to  find  a 
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real-rational  proper  K  to  minlalze  either  the  H2  or  the 
Hoo  norm  of  under  the  constraint  that  K  stabilizes  G.  T*<i 

is  the  transfer  function  from  e  to  d  in  Figure  2.4.  K 
stabilizes  G  if  the  states  of  G  and  K  tend  to  zero  for  every 
initial  condition  and  can  be  thought  of  as  a  simple 
definition  of  internal  stability.  The  H2  norm  is  given  by 
(2:831) 


(4r  J  t 
'■  -00 


.  l/'Z 


trlG(J«)  G(3M)))d« 


The  Hoo  norm  is  given  by  (2:631) 

loi.  = 

The  standard  problem  will  be  solved  using  the 
Doyle/Glover  parameterization(2) .  In  the  case  of  H2 
optimization,  the  following  is  required: 


(2.9) 


(2.10) 


1)  Dll  »  0  and  D22  »  0 

2)  D12’'*D12  »  I  and  021*021*^  =  I 

3)  (A,B1)  stabilizable  and  (Cl, A)  detectable 

4)  (A,B2)  stabilizable  and  (C2,A)  detectable 

To  use  the  Doyle/Glover  parameterization,  the  plant  input  and 
output  must  be  scaled  as  shown  in  Figure  2.5. 

The  following  equation  development  is  ta)cen  from  a 
summary  of  the  Doyle/Glover  parameterization( 2 )  by 
Ridgely(5:l-6) .  The  resulting  equations  from  the  scaling 
are 

u  -  8^'^  u  (2.11) 
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Fig*  2.5  Scaled  Standard  Block  Diagraa  for 
H2  Optiaization 
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(2.12) 


y  =  Sy  y 

The  scale  factors  are  found  using  a  Cholesky  decomposition 


of  the  right  hand  side  of  the  equations 

T  T 

Su^u  “  °12  ^12 


-1  -1  T 

Sy  (Sy  1)'' 


^21^21 


The  plant  equations  now  become 


Ax  +  Bj^d  +  B2U 


(2.13) 

(2.14) 

(2.15) 


where 


e  =  C^x  + 


y  =  C,x  +  Dj^d 


Bo  *  BoS 

2  2  U 


Co  =  S  Co 
2  y  2 


-1 


Di2  =  DjjSy 


^21  '  ®y^21 


-1 


(2.16) 

(2.17) 

(2.18) 

(2.19) 

(2.20) 
(2.21) 


The  family  of  all  H2  compensators  is  given  by 


•^ed  ' 

K  -  F^(J,Q) 
J 


-1. 


C^(sl  -  A^)  *Bj  +  Dj 


A.  «  A  -  K,Co  -  BoK^ 
j  f  2  2  C 


-  fK,  K,,I 


(2.22) 

(2.23) 

(2.24) 

(2.25) 

(2.26) 
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(2.27) 


Dj  -  [  J  J  ]  (2.28) 

•'e  '  ®2%  ♦  ^12^‘=1 

=  ■<®21®l'^^2  *  ^2^  (2.30) 

Kf  =  ^2^2*^  ■*■  (2.31) 


•'£1  =  ’'2'=l’'“  12  *  =2  *2- 

where  and  are  the  solutions  to  the  Ricattl  equations 

(A-B2Dj^2''^^1^'^^2*^2^^"®2^12'^S^‘^2®2®2^^2^^1^^1  *  °  (2.33) 

^^ifi  ^  ^  m 

<^-^21  ^2>^2""2<^-V21  -^2S%^2^®1®1  "  °  <2.34) 


where 


A 

<=1 

.  (I  -  Si25i2'')Ci 

(2.35) 

A 

®1 

.  Bjd  -  BjJSjj) 

(2.36) 

and  Q  is  any  stable,  strictly  proper,  transfer  function  such 


that 


|Q|J  <  -  ^ 


(2.37) 


is  the  optinal  2  norm  solution,  which  is  obtained 
when  Q  is  chosen  equal  to  zero.  A  subopt imal  solution  is 
obtained  for  Q  not  equal  to  zero. 

In  the  Doyle/Glover  paraaeter ization  for  Hao 
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1 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 


optimization,  the  following  conditions  must  be  satisfied: 

1)  Dll  =  0  and  D22  =  0 

2)  D12^*D12  =  I  and  D21*D2l'^  =  I 

3)  (A,B1)  controllable  and  (Cl, A)  observable 

4)  (A,B2)  stabilizable  and  (C2,A)  detectable 

In  addition,  the  parameterization  for  Hod  optimization  also 
requires  scaling  on  d  and  e  as  shown  in  Figure  2.6.  The  scale 
factors  are  given  by 


d 

=  d 

(2.38) 

e 

=  y  r  e 

(2.39) 

y 

=  S  y 

y  ^ 

(2.40) 

u 

-1  " 

=  "  U 

(2.41) 

where  again  Su  and  Sy  are  found  using  a  Cholesky 

decomposition  of  the  right  hand  side  of 

T 

S  S 

1  T 

“  -  D,.,  D, 

(2.42) 

u  u 

r  12  12 

and 

-  -1.  -l.T 

1  T 

(2.43) 

Sy  "(Sy  ^ 

*  ?  ^21^21 

The  plant  equations  become 

Si 

MM# 

■  Ax  +  Bj^d  +  B2U 

(2.44) 

e 

M  mf  00 

•  V  *  “ij” 

(2.45) 

y 

M  0*0  0*0 

-  CjX  +  Dj^d 

(2.46) 

where 
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®1 

=  /Ur 

®i 

(2.47) 

mr 

(2.48) 

^1 

=  /Ur 

^1 

(2.49) 

CM 

1  U 

-  iyCj 

(2.50) 

^12 

=  /Ur 

^12Su"' 

(2.51) 

*21 

=  /Ur 

"y''2i 

(2.52) 

The  family  of  all  stabilizing  compensators  such  that 


Is  given  by 


T 


K 


(2.53) 

M  A#  ^ 

ed 

F^(P,K) 

(2.54) 

K  = 

F^(J,Q) 

(2.55) 

J  = 

Cj(sl  -  Aj)“^B^  +  Dj 

(2.56) 

A  -  KjSj  -  ijKj  ♦ 

(2.57) 

=1  • 

(Kj  KjJ 

(2.58) 

■ 

[c] 

(2.59) 

■>j  - 

[?5] 

(2.60) 

K  - 
C 

~  T  ~  m~  _n 

(Bj  XOO  ♦  Dj2 

(2.61) 

cl 

-(D2iBi^X®  +  CjXI-YooX®)"^ 

(2.62) 

•'f  “ 

~  T  "  T 

YcCj^  .  BjDj/ 

(2.63) 
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(2.64) 


tl 


B. 


where  Xao  and  Yod  are  solutions  of  the  algebraic  Ricatti 
equations 

°  (2.65) 

and 


(A-Bj^D2j^’'c2)Y<»fYoo(A-Bj^D23^'‘’c2)’’-YocC2’’c2Yoi>fBj^Bj^^  =  0 

(2.66) 

where 

(2.67) 

Si  -  Bjd  -  £2i’'£2i) 

(2.66) 

and  Q  is  any  stable,  proper  transfer  function  such 

that 

|0|„  <  1 

(2.69) 

These  equations  are  solved  by  what  is  )cnovn  as  gamma 
iteration.  The  Hoo  optimal  solution  is  not  unique.  To  be 
valid,  Xoo  and  Yoo  must  have  solutions  that  are  elements  of  the 
Ricatti  domain  and  positive  semidef inite .  The  reasons  for 
this  are  covered  in  the  Doyle/Glover  development ( 2 : 833 ) .  The 
product  YooXoo  must  also  have  a  spectral  radius  less  than  one, 
where  the  spectral  radius  of  a  matrix  A  is  given  by 

P(A)  -  max|\, (A) I  (2.70) 

i  ‘ 

The  above  equations  were  solved  using  PC-HATLAB 
routines . 
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Ill .  PiQbien 


3.1  general  Description 

The  aim  of  the  design  was  to  develop  a  lateral  autopilot 
that  met  the  design  requirements.  The  autopilot  had  to  meet 
performance,  robustness,  control  activity,  and  gust  response 
requirements  over  a  range  of  fourteen  flight  conditions. 

The  autopilot  is  to  operate  in  two  modes:  heading  hold 
mode,  and  localizer  capture  and  track  mode.  In  the  heading 
control  mode,  the  autopilot  is  to  track  a  pilot-commanded 
heading  angle.  In  the  localizer  capture  and  track  mode,  it 
must  drive  the  cross  track  error  to  zero  and  align  the 
aircraft  velocity  vector  with  the  runway  heading.  Figure  3.1 
shows  a  diagram  of  the  approacn  geometry. 

In  this  thesis,  only  the  localizer  capture  and  track 
autopilot  will  be  designed.  The  localizer  capture  and  track 
autopilot  Is  more  difficult  to  design  since  It  must  track  two 
signals . 

3.2  Plant  Description 

The  state  space  representation  for  the  fourteen  flight 
conditions  was  provided.  These  are  in  the  form 

=  Ax  -t-  Bu  (3.1) 

y  *  Cx  +  Du  (3.2) 

The  state  space  representations  are  for  the  lateral  dynamics 
of  the  aircraft.  The  states  are 

*»(/?p^ryyl'^  (3.3) 
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Fig.  3.1  Approach  Oeoaetry 


where 

p  -  side  slip  angle  (deg) 
p  -  roll  rate  (deg/sec) 

-  roll  angle  (deg) 
r  -  yaw  rate  (deg/sec) 
y/  -  heading  angle  (deg) 
y  -  cross  track  error  (ft) 

The  outputs  are 

y  ~  I  p  lly  t  (3.4) 

where 

ft  -  side  slip  angular  rate  (deg/sec) 
n  -  side  acceleration  (g's) 

y  2 

t  -  yaw  acceleration  (deg/sec  ) 

-  ground  track  heading  angle 

The  control  Inputs  are 

u  =  I  6a  dr  j’’  (3.5) 

where 

6a  =  aileron  deflection  angle  (deg) 

6r  s  rudder  deflection  angle  (deg) 

The  disturbance  input  was 

pq  =  beta  gust  (deg) 

The  numbering  for  the  flight  conditions  is  not  sequential 
since  they  were  taken  from  a  larger  data  base  of  linearized 
flight  conditions.  The  flight  conditions  provided  are  001, 
004,  017,  020,  033,  036,  049,  052,  065,  068,  081,  084,  097, 
and  100.  The  state  space  representations  for  these  are  in 
Appendix  A. 
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3.3  Design  Requirements 


The  design  requirements  Include  specifications  on 
performance,  robustness,  control  activity,  and  gust  response. 

The  performance  specifications  required  zero  time  domain 
overshoot  and  zero  steady  state  tracking  error.  The  bank 
angle  could  not  exceed  30  degrees.  The  lateral  acceleration 
during  aircraft  maneuvers  was  to  remain  below  0.05  g*s.  The 
bandwidth  requirements  are  just  stated  as  to  maximize 
response  bandwldths  while  remaining  within  limitations  on 
overshoot  and  control  surface  activity. 

The  robustness  requirements  are  stated  In  terms  of  the 
eigenvalues  or  pole  placement.  A  minimum  damping  ratio  of 
0.4  Is  required  for  all  modes.  The  minimum  damping  ratio  on 
the  dominant  mode  Is  required  to  be  0.6.  Stability  margins 
of  4  db  and  40  degrees  simultaneously  are  required  for  all 
Inputs  and  outputs. 

The  control  surface  activity  for  the  ailerons  and 
rudder  must  remain  within  *  15  degrees.  The  control  surface 
rates  must  remain  below  +  30  deg/sec. 

There  were  additional  requirements  placed  on  the 
response  to  a  1  ft/sec  RMS  Dryden  gust.  The  heading  and 
sideslip  angular  positions  are  to  remain  below  0.5  degrees. 
The  angular  rates  for  these  two  states  must  remain  below  0.5 
deg/sec.  The  control  surface  activity  for  the  gust  response 
is  to  be  less  than  2  degrees  in  position  and  5  deg/sec  In 
rate  for  the  ailerons  and  rudder. 


19 


3.4  Open  Loop  Plant  Description 


The  flight  conditions  given  are  for  various  phases  of 
the  approach,  including  different  aircraft  configurations. 
The  data  can  be  broken  down  into  two  groups  based  on 
velocity.  Plight  Conditions  001,  017,  033,  049,  065,  081, 
and  097  are  for  a  slow  approach.  Plight  Conditions  004, 

020,  036,  052,  066,  084,  and  100  are  for  a  fast  approach. 
Table  I  list  the  aircraft  parameters. 

The  typical  lateral  dynamics  include  the  spiral,  roll, 
and  Dutch  roll  modes.  The  spiral  and  roll  modes  have  real 
roots  and  the  Dutch  roll  is  a  complex  conjugate  pair.  The 
roots  for  the  three  modes  for  Plight  Condition  084  are 

spiral:  =  -0.0107  (  Tj^  =  93.5  sec  ) 

roll:  *  -3.4824  (Tj  =  0.287  sec) 

Dutch  roll:  X_  .  =  -0.1278  +  1.5028 
3,4  • 

The  Dutch  roll  has  a  natural  frequency,  vn,  of  1.508  rad/sec 

and  is  very  lightly  damped  vith  a  damping  ratio  of  0.085. 

All  flight  conditions,  vith  the  exceptions  of  065  and  097, 
have  stable  roots.  The  characteristics  modes  are  given  in 
Table  II. 

The  heading  and  cross  track  error  states  were  added  for 
the  localizer  capture  and  track  mode.  These  two  states 
Introduced  tvo  poles  at  the  origin.  The  heading  angle  is 
given  by 


and  the  cross  track  by 
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Table  II  Characteristic 
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(3.7) 


y  =  ;  Vsin(v/  *  ft  -  ^ina) 

where 

V  -  aircraft  velocity  (ft/sec) 
o  -  angle  of  attack  (deg) 

An  explanation  of  the  localizer  capture  and  track  node  is 
given  in  Roskam  (6:1198-1207). 

To  evaluate  the  open-loop  frequency  domain 
character istics/  the  open-loop  singular  values  were  plotted. 
The  singular  value  plot  for  flight  conditions  001  and  004 
are  shown  in  Figures  3.2  -  3.3.  The  general  trends 
are  that  the  plant  has  widely  separated  singular  values  at 
low  frequency.  The  plant  doesn't  exhibit  integral  action. 
The  a  has  very  low  gain  at  low  frequency.  This  will  need  to 
be  corrected  for  good  performance  and  steady  state  tracking 
accuracies.  The  plant/  as  given/  is  proper  and  doesn't  roll 
off  at  high  frequency.  This  will  need  to  be  changed  for 
robustness/  and  to  be  able  to  use  H2  and  Hop  optimization. 
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FREQ  -  rad/s 

Fig.  3.2  Flight  Condition  001  Open  Loop  Singular  Value  Plot 
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IV.  PiQbleffi  Synthesis 

4.1  Modification  ths.  Plant 

The  state  space  description  of  the  plant  as  given  is 
not  strictly  proper.  One  of  the  requirenents  of  H2  and  Hao 
theory  is  that  the  D22  matrix  of  equation  (2.8)  be  zero. 

This  means  that  the  plant  can  have  no  feedforward  elements. 

In  order  to  satisfy  this  requirement,  actuator  dynamics  were 
augmented  into  the  plant.  The  actuators  for  the  aileron  and 
rudder  were  represented  by 

6a (s)  =  6r(s)  =  (4.1) 

The  state  space  representation  of  the  actuators  is  given  by 

6  «=  X6  +  B6  (4.2) 

and 

y  =  C6  (4.3) 

where 

[-i-J]  ■>-  [S  1]  c.  [^5  J] 

The  actuator  dynamics  were  augmented  into  the  plant  by 
cascading  the  two  systems  together  as  shown  in  Figure  4.1. 


Fig.  4.1  Cascaded  Systems 

This  yields  a  new  state  space  representation  with  eight 
states.  The  singular  values  were  replotted  for  flight 
conditions  001  and  004,  and  as  expected  for  a  strictly  proper 
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systen,  the  singular  values  rolled  o££  at  20  db/decade  at 
high  £requency  as  shovn  in  Figures  4.2  and  4.3. 

In  addition  to  augmenting  the  actuator  dynamics,  the 
cross  track  error  was  scaled  by  100  to  normalize  outputs. 

Nov  the  system  could  be  £ormulated  Into  the  standard 
problem  shovn  In  Figure  2.4. 

4.2  gmaJLl  Gain  Problem  Formulation 

In  order  to  £ormulate  the  control  design  Into  the 
standard  problem,  a  suitable  block  diagram  vas  needed.  The 
block  diagram  In  Figure  4.4  vas  used. 

The  system  has  exogenous  Inputs,  d,  o£  tracking  signals, 
gust  Input,  and  sensor  noise.  These  are  labeled  dl,  d2,  and 
d3  respectively.  The  tracking  signals  £or  the  localizer 
capture  and  hold  autopilot  are  the  heading  angle  and  cross 
track  error.  The  only  disturbance  Is  the  gust  Input. 

The  outputs  o£  the  system  as  given  vere 

V  .  t  P  nj,  t  (4.4) 

These  vere  mod 1£ led  to 

V  =  t  ^  "y  v  y  j’’ 

The  output  of  the  cross  track  error  vas  added  since  the  error 
signal  for  this  state  Is  desired.  The  yav  acceleration  vas 
deleted  since  It  vas  not  needed. 

The  dl  vector  vas  made  the  same  as  v  to  keep  the  number 
of  Inputs  and  outputs  equal.  This  vas  done  In  order  to 
formulate  the  problem  vlth  a  tracking  error  signal  of  r-v  as 
is  done  in  the  tracking  problem  setup  in  Francis  (3:3. 4-3. 5). 
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Fig.  4.4  Systen  Block  Diagram 


(4.6) 


The  control  input  Is  also  included  in  the  cost  function 

[l  r-v|^  ♦ 

where  p  is  a  weighting  factor  to  penalize  the  amount  of 
control  power  used. 

Weights  were  put  on  the  incoming  signal,  dl,  and  the 
error  signal,  el.  These  were  added  to  make  the  system 
tunable  when  designing  the  compensators  and  to  satisfy  the 
H2  and  Hoo  requirements. 

The  system  block  diagram  in  Figure  4.4  is  a  hybrid  of 
the  plant  state  space  representation  and  the  transfer 
functions  of  the  weights  that  have  been  added.  To  put  the 
system  in  the  form  of  equations  (2.6)-(2.6),  the  weights  need 
to  be  written  in  a  state  space  form.  These,  along  with  the 
plant  state  space  representation,  are  then  substituted  into 
equations  (2.6)-(2.8). 


The  weight  on  dl,  Wl,  can  be  written  as 

^1  ~  ^1^1  ^  ®1^1 

(4.7) 

r  . 

(4.8) 

The 

weight  on  the  error.  We,  can  be  written  as 

(4.9) 

®1  '  Ve  * 

(4.10) 

The 

weight  on  the  control  input  is  a  diagonal  scale 

factor 

given  by 

«2  ‘  ^ 

Mow  using  Figure  4.2, 

(4.11) 

P  P  P  g  2  p 

(4.12) 
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^  P  P 

(4.13) 

''  =  S^p  "  V2  "  ""3^3 

(4.14) 

These  equations  can 

be  used  to  write 

*e  *  *e*e*  * 

B  D-d,  -  BCx  -  BDd.,-B  D-d- 
ell  epp  eg2  ejj 

(4.15) 

*1  “  * 

D^D,d-  -  D^C  -  D^D  d-  -  D^D-d- 

eii  epp  eg2  ejj 

(4.16) 

The  standard  problem  state  space  form  can  now  be 
written  for  Figure  4.4  as 
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L«3J 


[u] 

(4.17) 

(u) 

(4.18) 


+  [0  ]lul 

(4.19) 


The  conditions  for  the  Doyle/Glover  parameterization 
require  Dll  and  D22  to  be  zero.  D22  has  already  been  ta)cen 
care  of.  Dll  is  made  equal  to  zero  by  choosing  the  weights 
properly.  From  equation  (4.16),  it  can  be  seen  that  D^  must 
be  zero.  This  says  that  the  weight  on  the  error  signal  must 
be  strictly  proper. 

The  other  conditions  are  that  D12^D12  and  D21D21^  be 
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similar  to  I.  The  first  condition  Is  satisfied.  In  that 

D12'^D12  *  (4.20) 

which  Is  chosen  to  be  a  full  ran)c  diagonal  matrix.  The  later 
Is  satisfied  by 

D21D2l'’^  =  +  D-D-’’  (4.21) 

11  9  9  J  J 

where  Is  a  4  by  4  Identity  matrix  added  for  Inputlng 
measurement  noise,  Is  a  4  by  1  vector  the  gust 
disturbance  enters  the  output  of  the  system  through,  and  D1 
is  a  4  by  4  matrix  since  W1  is  proper.  This  guarantees  that 
equations  (4.20)  and  (4.21)  have  full  ran)c. 

The  open  loop  plant  Is  both  controllable  and  observable. 
The  augmented  system  Is  stabilizable  and  detectable. 

Condition  3  for  Hoo  requires  that  (A,B1)  be  controllable  and 
(A, Cl)  be  observable.  These  conditions  are  not  truly  met,  as 
the  condition  numbers  of  the  controllability  and 
observability  matrices  are  very  large.  The  condition  number 
Is  the  ratio  of  the  maximum  and  minimum  singular  values.  A 
large  condition  number  Is  Indicative  of  a  nearly  singular 
matrix.  This  will  probably  cause  numerical  problems  In  the 
computation  of  the  Hoo  controllers.  This  problem  will  be 
addressed  later  In  the  discussion  section. 

The  system  equations  can  be  modified  now  that  the  form 
of  the  weights  has  been  specified.  The  new  equations  are 
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(u] 

(4.23) 


[y]“C-S 


+  [0  ]lu) 

(4.24) 


The  problem  synthesis  is  complete.  These  equations  vill 
nov  be  used  in  the  H2  and  Hoo  routines  to  design  two 
controllers  for  the  aircraft's  lateral  autopilot. 
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V.  Analysis 


5.1  General  Approach 

The  goal  of  this  thesis  was  to  design  one  controller 
that  Beets  all  the  design  requireaents  for  all  flight 
conditions.  It  was  determined  that  one  controller  would  not 
meet  the  specifications  for  all  flight  conditions  after  a 
controller  designed  for  a  nominal  case  for  all  fourteen 
flight  conditions  yielded  poor  results  in  closed  loop 
simulation  runs.  The  flight  conditions  were  then  broken  down 
into  the  slow  and  fast  conditions.  A  nominal  condition  was 
picked  for  which  the  design  was  based  on.  Once  these 
controllers  were  designed,  they  were  used  in  the  closed  loop 
simulation  to  evaluate  the  performance.  A  controller  was 
obtained  using  H2  and  Hoo  control  theory  for  both  the  slow  and 
fast  approach  flight  conditions. 

The  techniques  used  to  design  each  controller  was  a  mix 
of  trying  to  shape  the  singular  value  plots  of  the  loop  gain, 
the  sensitivity  function,  and  the  complementary  sensitivity 
function  combined  with  evaluating  the  time  domain  response 
characteristics.  The  weights  on  the  input  signal,  Wl,  and  on 
the  error  signal.  We,  were  varied  in  trying  to  obtain  the 
desired  response.  In  addition  to  these,  the  scale  parameter, 
Pf  on  the  control  Inputs  was  also  varied. 

The  program  PCMATLAB  was  used  for  this  thesis.  Many 
programs  were  written  using  the  PCMATLAB-based  software.  The 
Control  System  Toolbox  was  also  used. 
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Routines  to  build  the  natrices  in  equations  (4.22)> 
(4.24)  vere  written.  DOABCD  formulated  the  A,  Bl,  B2,  Cl, 

C2,  Dll,  D12,  D21,  and  the  D22  matrices  for  input  into  the 
code  to  solve  for  the  H2  and  the  Hoo  controllers.  These 
programs  returned  the  H2  and  Hod  controllers . 

Routines  to  build  the  state  space  matrices  for  the 
weights  were  also  written.  These  routines,  DOWl  and  DOWE, 
returned  the  Al,  Bl,  Cl,  Dl,  Ae,  Be,  and  Ce  matrices  used  in 
equations  (4.22)-(4.24) . 

The  routine  CLOSE  formed  the  closed  loop  matrices 
which  consisted  of  the  original  plant  and  the  controller. 

The  PCMATLAB  routine  LSIM  was  used  to  simulate  the  closed 
loop  system. 

Other  utility  routines  were  written.  The  routines  RUNH2 
and  RUNHI  used  the  above  mentioned  routines  to  do  the  design 
form  start  to  finish.  Several  routines  to  plot  simulation 
data  and  singular  value  plots  vere  also  written.  The  user 
written  code  is  listed  in  Appendix  B. 

5.2  Nominal  Controller  QfiAign  And.  AnalYSig 

The  flight  condition  chosen  for  the  nominal  controller 
design  for  the  fast  approach  data  was  Flight  Condition  084. 
Flight  Condition  033  was  chosen  for  the  nominal  controller 
design  for  the  slow  approach  data.  The  data  is  a  sampling  of 
the  most  extreme  conditions  for  the  landing  profile. 

The  weights  used  in  the  design  vere  based  on  the 
classical  control  ideas  of  bandvidths  and  gains.  The 
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response  could  be  altered  significantly  vith  the  weights. 

Weighting  strategies  of  trying  to  shape  the  singular 
value  plots  were  also  used.  This  entailed  knowing  what  the 
sensitivity  of  the  system  should  look  like  and  using  a  weight 
which  gave  the  desired  sensitivity.  The  sensitivity  should 
be  low  at  low  frequency  since  the  loop  gain  is  large.  The 
sensitivity  function  is  given  by 

S  =  (I  +  GK)“^  (5.1) 

On  a  Bode  plot,  the  sensitivity  goes  to  0  db  at  high 
frequency  since  GK  is  small.  The  weights  were  also  chosen  in 
an  effort  to  achieve  high  loop  gains  at  low  frequency. 

The  Bode  magnitude  plot  of  Wl,  the  weight  on  the  inputs, 
is  shown  in  Figure  5.1.  This  weight  is  on  the  ground 
track  heading  angle  and  the  cross  track  error.  The  other 
inputs  were  just  passed  through  with  unity  gain  since  they 
were  mainly  there  to  achieve  an  equal  number  of  inputs  to  the 
number  of  outputs. 

The  W1  weights  were  chosen  to  shape  the  inputs  of  the 
system  as  a  step.  Moving  the  bandwidths  had  significant 
effect  on  the  system  response,  as  did  changing  the  DC  gains. 

The  error  weights.  We,  were  just  as  effective.  The  error 
weights  were  thought  of  in  terms  of  penalizing  a  particular 
signal.  A  good  example  of  this  is  the  lightly  damped  pole 
from  the  Dutch  roll  mode  around  1.5  rad/sec.  This  pole  had  a 
damping  ratio  of  0.08,  for  Flight  Condition  064,  and  was 
increased  to  0.52  by  penalizing  the  sideslip  angular  rate 
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error,  The  Bode  plots  of  the  weights  for  the  four  error 

outputs  are  shown  in  Figures  5.2-5. 5.  The  weights  on  ft  and 
n^  each  added  one  state  to  the  compensator.  The  choice  of 
weights  on  the  ground  track  heading  error,  r  and  cross 

track  error,  y^,  were  made  in  terms  of  adjusting  the 
gains  for  steady  state  accuracy  and  moving  the  bandwidths  for 
damping  the  response. 

The  flexibility  of  the  weights  gave  enormous  amounts  of 
tuning  control.  Unfortunately,  it  also  gives  an  enormous 
amount  of  possible  combinations  of  weighting.  As  mentioned, 
several  weighting  strategies  were  used.  The  resultant 
strategy  was  somewhat  of  a  combination  of  all  of  them. 

The  success  of  shaping  the  singular  value  loop  shapes 
was  minimal.  The  minimum  singular  value  of  the  loop  gain  was 
very  low.  This  caused  problems  in  obtaining  the  desired 
steady  state  accuracy  in  both  channels. 

When  the  weights  were  finalized  for  the  nominal 
controllers,  the  response  for  the  closed  loop  system  was  run 
for  all  flight  conditions.  This  will  be  discussed  in  the  next 
section.  The  closed  loop  system  was  formed  by  writing  the 
controller  in  its  state  space  form,  given  by 

*k  *  ^k’^k  *  (5.2) 

u  «  Cj^Xj^  (5.3) 

Closing  the  loop  yielded  the  equations 


39 


100 


40 


Weight  on  Betadot  Bcror  Output 
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Pig.  5.3  Weight  on  Noraal  Acceleration  Error  Output 
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Fig.  5.4  Weight  on  Ground  Track  Heading  Angle  Error  Output 
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Fig.  5.S  Weight  on  Cross  Track  Error  Error  Output 
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The  controller  returned  by  the  H2  and  Hoo  calculations  has 
the  number  o£  states  equal  to  the  original  plant  plus  the 
number  of  the  states  associated  vith  the  weighting  filters. 
The  final  controller  has  16  states. 

The  initial  conditions  input  to  the  simulation  were  an 
initial  cross  track  error  offset  of  2500  feet.  This 
corresponded  to  intercepting  a  2.5  degree  localizer  beam 
at  9.5  nautical  miles  from  the  transmitter.  The  initial 
heading  angle  was  aligned  vith  the  runway  heading  angle  which 
was  0  degrees. 

The  response  obtained  for  the  nominal  system  for  the 
H2  controller  designed  using  Plight  Condition  084  and  033  are 
shown  in  Figures  5.6  and  5.7.  The  response  for  these 
controllers  met  all  the  specifications.  All  poles  had  a 
damping  ratio  greater  than  the  required  0.4  and  the  dominant 
pole  was  greater  than  0.6.  The  steady  state  tracking  for  the 
cross  track  error  and  ground  track  heading  angle  also  met  the 
specifications.  The  control  surface  angular  positions  and 
rates  were  within  the  allowable  values. 
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Plq.  5.7a.-  d.  Pit  Cond  033  Tlae  Response  for  H2  Ctrl 
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Fig.  5.7e.-  £.  Pit  Cond  033  Tl*e  Response  for  H2  Ctrl 


The  singular  value  plots  for  Flight  Condition  084  and 
033  are  shown  in  Figures  5.8  -  5.9.  The  loop  gain  for  both 
cases  is  not  high  enough  for  good  performance  robustness. 

This  will  be  addressed  in  the  next  section  when  these  two 
controllers  are  used  for  the  off  nominal  conditions.  The 
maximum  singular  value  of  the  closed  loop  or  complementary 
sensitive  function  for  both  cases  indicates  adequate  stability 
margins . 

The  approach  for  the  Hoo  controller  design  was  basically 
the  same  as  for  the  H2  controller  except  for  the  addition  of 
the  gamma  iteration.  Gamma,  the  scale  factor  shown  in 
Figure  2.6,  must  be  chosen  close  to  the  optimum  Hoo  norm. 

This  is  necessary  since  the  Hoo  optimal  controller  can  only  be 
approached  and  not  calculated  directly  as  in  H2  optimization. 
The  weights  used  for  the  H2  controller  design  required 
modifications  to  allow  the  Hoo  computations  to  achieve 
reasonable  results.  By  this,  it  is  meant  that  initially  the 
controller  developed  using  Hoo  control  theory  and  Flight 
Condition  084  destabilized  the  closed  loop  system.  The 
reason  for  the  destabilizing  controllers  was  numerical  error. 
The  controllability  and  observability  conditions  required  for 
(A,B1)  and  (Cl, A)  in  Hoo  control  theory  were  not  met. 

The  initial  value  of  gamma  used  to  achieve  a  solution 
was  170.  This  was  improved  upon  by  bringing  the  separation 
of  the  poles  and  zeros  in  the  weights  closer  together.  The 
final  controller  designed  for  the  fast  approach  data  using 
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Gain  SV  .  b.  SenaiiiTity  Func 


Flq.  5.8a.-  c.  Fit  Cond  084  Singular  Value  Plots  for  H2  Ctrl 


Flight  Condition  084,  had  a  gamma  of  67  and  a  spectral  radius 
of  0.36.  Similar  problems  were  encountered  for  Flight 
Condition  033.  The  controller  designed  for  this  case  had  a 
gamma  and  rho  of  77  and  0.47,  respectively.  The 
modified  weights  allowed  gamma  to  be  reduced  and  the  optimum 
Hod  norm  to  be  approached,  but  not  close  enough  to  achieve 
desirable  results.  The  maximum  Hoo  norm  for  T^^  was 
calculated  to  be  about  0.85.  More  discussion  about  the 
reasons  for  some  of  the  results  obtained  will  follow  in  the 
next  section. 

The  final  designs  using  Hoo  optimization  yielded  results 
that  met  the  specifications  for  the  nominal  flight  conditions. 
The  time  plots  for  the  two  cases  are  shown  in  Figures  5.10 
and  5.11.  The  tracking  and  response  of  the  system  was 
acceptable.  The  control  surface  activity  was  within  the 
specified  limits.  The  damping  for  all  modes  was  greater  than 
0.4  and  the  dominant  mode  damping  was  greater  than  0.6.  The 
roll  angles  and  normal  acceleration  were  well  within  the 
limits. 

The  singular  value  plots  for  the  Hoo  designs  for  Flight 
Condition  084  and  Flight  Condition  033  are  shown  in  Figures 
5.12  and  5.13.  The  loop  gain  minimum  singular  values  show 
that  the  system  does  not  have  adequate  gain  at  low 
frequency  which  implies  the  system  will  not  be  robust  to 
parameter  variations.  The  bandwidth  of  the  system  is  also 
very  low. 
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Response  for  Ha>  Ctrl 
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Flq.  5.13a.-  c.  Pit  Cond  033  Singular  Value  Plots  for  Ho»  Ctrl 


VI .  Results 


6.1  Closed  Loop  Simulation  Results 

The  closed  loop  simulation  was  run  using  the  controllers 
developed  using  the  H2  and  Hoo  control  theory.  The  H2 
compensator  developed  for  Flight  Condition  084  was  used  as 
the  controller  for  flight  conditions  004,  020,  036,  052,  068, 
084,  and  100,  the  fast  approach  flight  conditions.  The  H2 
compensator  designed  for  Flight  Condition  033  was  used  in  a 
closed  loop  simulation  for  flight  conditions  001,  017,  033, 
049,  065,  081,  and  097,  the  slow  approach  flight  conditions. 
The  time  plots  for  the  H2  simulation  results  are  shown  in 
Figures  6.1-6.12. 

The  simulation  results  for  the  off  nominal  conditions 
were  varied.  Flight  Conditions  001  and  004  had  an  unstable 
pole.  Looking  at  the  time  history  plots  for  these  two  cases, 
it  is  clearly  evident  that  Flight  Condition  001  is  unstable. 
The  same  is  not  true  for  Flight  Condition  004.  The  time 
constant  for  the  unstable  mode  is  very  large  and  does  not 
show  up  over  the  duration  of  time  plotted.  All  other  flight 
conditions  were  stable  and  had  a  damping  ratio  of  greater 
than  0.6  for  the  dominant  mode.  Flight  Conditions  065 
through  100,  with  the  exception  of  084,  had  modes  with 
a  damping  ratio  of  less  than  0.4.  The  lowest  damped  mode  for 
these  cases  was  0.27. 

The  steady  state  tracking  accuracies  were  also  varied. 
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Fig.  6.3a.-  d.  Fit  Cond  036  Tlae  Response  for  H2  Ctrl 
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Pig.  6.6a.-  d.  Pit  Cond  100  Tlae  Response  for  H2  Ctrl 
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Pig.  6.8e.-  £.  Pit  Cond  017  Tlae  Response  for  H2  Ctrl 
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The  fast  approach  flight  conditions  tracked  reasonably  veil. 
Several  of  the  cases  had  not  achieved  zero  steady  state 
tracking  accuracies  during  the  desired  time  interval.  The 
slow  approach  cases  took  even  longer  in  approaching  the 
desired  steady  state  values. 

The  control  surface  activity  was  well  within  the 
specifications.  The  angular  positions  and  rates  were  within 
±15  degrees  and  the  ±30  deg/sec,  respectively.  The  bank 
angle  never  exceeded  about  ±17  degrees.  The  lateral 
acceleration  never  exceeded  0.01  g's. 

The  Hoo  results  were  not  as  favorable  when  the 
compensators  designed  for  Flight  Conditions  033  and  084  were 
used  in  the  closed  loop  simulation  for  the  other  flight 
conditions.  The  results  for  the  fast  approach  flight 
conditions  yielded  an  unstable  mode  in  all  cases  except  for 
Flight  condition  004.  The  results  for  the  slow  approach  data 
had  an  unstable  mode  in  all  cases. 

6.2  Discussion  fii.  ReswUs 

The  performance  for  the  controllers  developed  did 
not  achieve  the  performance  desired.  Originally,  the  goal 
was  to  design  one  controller  using  H2  and  Hoo  control  theory 
that  would  meet  the  performance,  stability,  and  robustness 
requirements.  It  was  found  that  the  velocity  dependence  of 
the  performance  required  the  data  to  be  grouped  into  two  sets 
of  data.  It  was  hoped  that  controllers  could  then  be 
developed  that  would  yield  acceptable  results.  The  intent 
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was  then  to  compare  the  results  of  the  H2  and  Hoo  results. 

The  H2  and  Hod  control  theory  did  yield  acceptable 
results  for  the  nominal  cases  for  which  the  compensators  were 
designed.  All  specifications  were  met  for  the  nominal  fast 
and  slow  approach  flight  conditions. 

The  H2  compensators  did  yield  better  performance  when 
used  in  the  off  nominal  cases.  This  was  expected  since  the 
loop  gain  for  the  H2  designs  was  higher  therefore  implying 
better  performance.  Adequate  loop  gains  still  could  not  be 
achieved  for  either  compensator  with  H2  or  Hoo  control  theory. 

It  was  noted  earlier  that  the  controllability  and 
observability  matrices  were  poorly  conditioned.  The 
requirements  for  H2  and  Hoo  control  theory  are  that  the  (A,B1) 
and  the  (A, Cl)  matrices  be  stabilizable  and  detectable  for  H2 
and  controllable  and  observable  for  Hoo.  The  poor 
conditioning  of  these  matrices  obviously  resulted  in  numerical 
problems,  since  these  requirements  were  violated. 

The  effects  were  clearly  noticeable  when  designing  the 
Hoo  controllers.  The  weights  used  for  the  H2  design  did  not 
work  for  the  Hoo  design.  The  value  of  gamma  required  in  the 
gamma  iteration  yielded  an  Hoo  controller  that  was  very  close 
to  the  H2  controller.  This  occurred  because  as  gamma 
approaches  a  large  number,  the  Hoo  equations  reduce  to  the  H2 
equations.  There  was  also  some  noticable  effects  when  the 
weights  were  changed.  As  the  weights  were  changed,  checking 
the  controllability  and  observability  matrices  revealed  that 
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the  conditioning,  while  still  very  poor,  was  Improving.  The 
condition  number  was  reduced  by  10  orders  of  magnitude  by 
adjusting  the  weights.  As  mentioned  before,  the  results  were 
that  the  value  of  gamma  could  be  reduced  and  the  spectral 
radius  Increased. 

The  poor  conditioning  of  the  plant  may  have  been 
responsible  for  not  being  able  to  achieve  good  control 
response.  An  analysis  of  the  plant  conditioning  revealed  It 
was  poor.  In  (1:586),  it  Is  stated  that  "  the  plant  condition 
number  provides  a  quantitative  means  for  assessing  the  degree 
of  difficulty  In  achieving  good  control.”  In  addition,  the 
comment  is  made  that  to  obtain  a  good  loop  shape  (high  gain 
at  low  frequency),  while  ensuring  robust  stability,  the  plant 
must  be  well  conditioned. 

The  selection  of  weights  was  very  difficult.  The 
process  was  highly  Iterative.  The  process  of  selecting 
weights,  checltlng  the  singular  value  plots,  and  loo)(lng  at 
the  time  response  plots  was  very  time  consuming.  The  weights 
selected  affected  the  response  dramatically.  Even  with  all 
this  ability  to  shape  the  response,  the  final  selection  of 
weights  did  not  yield  the  desired  performance. 
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VII.  Conclusions  aM  Recommendations 


7.1  Summary  and  Conclusions 

The  objective  of  this  thesis  was  to  design  a  lateral 
autopilot  for  a  transport  aircraft  using  H2  and  Hoo  control 
theory.  The  autopilot  was  to  be  robust  with  respect  to 
variations  In  flight  conditions. 

After  background  Information  was  covered  In  Chapter  1, 
the  H2  and  Hod  control  theory  was  covered.  Important  aspects 
of  multivariable  control  theory  were  also  discussed.  The 
standard  problem  formulation  vas  Identified. 

In  Chapter  3  the  problem  vas  stated.  The  data 
provided  Included  the  state  space  representations  of  fourteen 
flight  conditions  In  addition  to  specifications  on 
performance,  robustness,  and  control  activity. 

This  data  vas  then  taken  and  formulated  Into  the 
standard  problem  or  the  small  gain  problem.  Once  In  this 
form,  an  Iterative  method  of  design  techniques  vas  used  in 
developing  a  controllers  using  H2  and  Hod  theory. 

In  Chapter  5,  the  controller  development  vas  covered  and 
the  performance  for  the  nominal  cases  vas  given.  Then  In 
Chapter  6,  the  controllers  developed  for  the  nominal  cases 
vere  used  in  a  closed  loop  simulation  and  evaluated  over 
fourteen  flight  conditions. 

The  controllers  developed  using  the  H2  control  theory 
provided  the  best  performance  and  stability  robustness  for 
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the  flight  conditions  tested.  While  these  controllers 
provided  some  measure  of  stability  and  performance 
robustness,  they  still  did  not  achieve  the  desired  results. 
The  111  conditioning  of  the  plant  may  have  limited  the 
achievable  results. 

The  Hcd  controllers  did  not  stabilize  the  off  nominal 
conditions.  The  (A,B1)  and  (A, Cl)  matrices  being 
uncontrollable  and  unobservable  numerically  Is  the  apparent 
reason  for  the  poor  performance  of  the  Hod  controllers.  This 
led  to  numerical  problems  as  evidenced  by  the  destabilizing 
controllers.  The  optimal  Hoo  solution  could  not  be 
approached.  The  minimum  value  of  gamma  that  could  be  used 
was  veil  above  the  optimum  value. 

The  controllers  obtained  using  Hoo  control  theory  were 
different  from  the  H2  controllers,  but  they  were  not 
optimum.  Due  to  the  numerical  errors  In  the  calculations. 

It  Is  unclear  exactly  hov  to  characterize  them. 

The  velght  selection  provided  a  large  amount  of 
flexibility  In  shaping  the  system  response,  but  the  desired 
system  response  and  robustness  vere  not  achieved.  An 
extensive  Iterative  process  of  velght  selection  vas 
accomplished.  The  results  given  vere  the  best  obtained.  It 
Is  thought  that,  because  of  the  extensive  search,  the 
required  performance  may  not  be  achievable  from  velght 
selection  alone  for  this  problem  setup. 
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7.2  Recommendations 


The  problem  synthesis  and  weight  selection  consumed  the 
majority  of  the  effort.  The  process  of  setting  up  the 
problem  is  very  undefined  but  is  also  very  important.  This 
problem  as  given  had  a  poorly  conditioned  plant.  The  outputs 
chosen  resulted  in  very  low  loop  gains  for  the  open  loop 
system.  A  different  set  of  outputs  can  yield  a  plant 
that  is  better  conditioned.  A  look  at  using  state  feedback 
shoved  a  considerable  improvement  in  the  open  loop  system. 

This  approach  was  not  used  since  ft  could  not  be  measured. 

If  the  plant  conditioning  can  be  eliminated  as  a 
concern/  the  numerical  problems  encountered  while  using  Hao 
control  theory  may  disappear  or  improve.  If  the  numerical 
difficulty  associated  with  Hoo  control  theory  is  still  present/ 
it  would  be  beneficial  to  quantify  what  the  numerical  limits 
on  the  conditioning  of  the  controllability  and  observability 
matrices  are  that  define!^  when  a  system  is  uncontrollable  and 
unobservable  numerically. 

This  leaves  the  area  of  weight  selection.  The  selection 
of  a  weighting  strategy  is  very  Important.  This  is  closely 
tied  in  with  the  problem  setup.  The  problv^iJ  t^s  setup  had 
weighting  filters  on  the  inputs  and  error  si^^  als.  As 
mentioned  previously/  the  weights  also  contributed  to  the 
numerical  problems.  In  trying  to  achieve  the  desired  loop 
shapes/  the  weight  selection  caused  the  conditioning  of  the 
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controllability  and  observability  matrices  of  (A,B1)  and 
(A, Cl)  to  be  poor,  respectively. 

The  process  of  weight  selection  is  highly  iterative  and 
time  consuming.  A  more  structured  way  of  choosing  weights  is 


needed . 


Appendix  A:  Lateral  sca£fi.  Data 

The  following  data  is  the  state  space  models  for  the 

fourteen  flight  conditions  for  the  Boeing  737  aircraft.  The 

data  is  presented  as  it  was  given.  The  states  are 

T 

x=t/3p0rv/y] 

The  outputs  are 

•  « 

y  -  (  Hy  r  ) 

The  inputs  to  the  system  are 

u  =  (  6a  6^  ) 

The  flight  conditions  are  001,  004,  017,  020,  033,  036,  049, 
052,  065,  066,  081,  064,  097,  and  100.  The  A,B  C,  and  D 
matrices  are  listed  by  letter  and  flight  condition  number. 


A-1 


AOOl 

-0.1481 

0.1984 

0.1355 

-0.9699 

0 

0 

-5.3738 

-2.0194 

-0.0025 

0.6359 

0 

0 

0 

1.0000 

0 

0.1990 

0 

0 

0.7480 

-0.1671 

-0.0023 

-0.1937 

0 

0 

0 

0 

0 

1.0196 

-0.0001 

0 

4.0649 

0 

-0.7932 

0 

4.0649 

-0.0002 

BOOl 

0.0019 

0.0474 

0.1481 

0.8651 

0.9043 

5.3738 

0 

0 

0 

0.0531 

-0.7397 

-0.7480 

0 

0 

0 

0 

0 

0 

cool 

-0.1481 

0.1984 

0.1355 

-0.9699 

0 

0 

-0.0187 

0.0004 

0 

0.0014 

0 

0 

0.7480 

-0.1671 

-0.0023 

-0.1937 

0 

0 

1.0000 

0 

-0.1951 

0 

1.0000 

0 

DOOl 

0.0019 

0.0474 

0.1481 

0.0002 

0.0060 

0.0187 

0.0531 

-0.7397 

-0.7480 

0 

0 

0 

A004 

-0.2130 

0.0645 

0.0906 

-0.9866 

0 

0 

-10.5010 

-3.0712 

-0.0086 

0.6389 

0 

0 

0 

1.0000 

0 

0.0654 

0 

0 

2.0574 

-0.1849 

-0.0031 

-0.2825 

0 

0 

0 

0 

0 

1.0021 

-0.0001 

0 

6.1857 

0 

-0.4035 

0 

6.1857 

-0.0002 

B004 

0.0035 

0.0686 

0.2130 

2.3464 

2.2789 

10.5010 

0 

0 

0 

0.1241 

-1.6316 

-2.0574 

0 

0 

0 

0 

0 

0 

C004 

-0.2130 

0.0645 

0.0906 

-0.9866 

0 

0 

-0.0410 

-0.0001 

0 

0.0022 

0 

0 

2.0574 

-0.1849 

-0.0031 

-0.2825 

0 

0 

1.0000 

0 

-0.0652 

0 

1.0000 

0 

D004 

0.0035 

0.0686 

0.2130 

0.0007 

0.0132 

0.0410 

0.1241 

-1.6316 

-2.0574 

0 

0 

0 

A-2 


I 


A017 

-0.1392 

0.1894 

0.1440 

-0.9714 

0 

0 

-4.9704 

-1.9664 

-0.0020 

0.7383 

0 

0 

0 

1.0000 

0 

0.1907 

0 

0 

0.6439 

-0.1689 

-0.0026 

-0.1771 

0 

0 

0 

0 

0 

1.0180 

-0.0001 

0 

3.8291 

0 

-0.7171 

0 

3.8291 

-0.0002 

B017 

0.0010 

0.0448 

0.1392 

0.9928 

0.8127 

4.9704 

0 

0 

0 

0.0843 

-0.6606 

-0.6439 

0 

0 

0 

0 

0 

0 

C017 

-0.1392 

0.1894 

0.1440 

-0.9714 

0 

0 

-0.0166 

0.0003 

0 

0.0013 

0 

0 

0.6439 

-0.1689 

-0.0026 

-0.1771 

0 

0 

1.0000 

0 

-0.1873 

0 

1.0000 

0 

D017 

0.0010 

0.0448 

0.1392 

0.0001 

0.0053 

0.0166 

0.0843 

-0.6606 

-0.6439 

0 

0 

0 

A020 

-0.2135 

0.0398 

0.0907 

-0.9878 

0 

0 

-10.4604 

-3.2899 

-0.0088 

0.7375 

0 

0 

0 

1.0000 

0 

0.0423 

0 

0 

2.1236 

-0.1767 

-0.0042 

-0.2651 

0 

0 

0 

0 

0 

1.0009 

-0.0001 

0 

6.1857 

0 

-0.2615 

0 

6.1857 

-0.0002 

B020 

0.0024 

0.0686 

0.2135 

2.3831 

2.3279 

10.4604 

0 

0 

0 

0.1444 

-1.6426 

-2.1236 

0 

0 

0 

0 

0 

0 

C020 

-0.2135 

0.0398 

0.0907 

-0.9878 

0 

0 

-0.0410 

-0.0005 

0 

0.0022 

0 

0 

2.1236 

-0.1767 

-0.0042 

-0.2651 

0 

0 

1.0000 

0 

-0.0423 

0 

1.0000 

0 

D020 

0.0024 

0.0686 

0.2135 

0.0005 

0.0132 

0.0410 

0.1444 

-1.6426 

-2.1236 

0 

0 

0 

A-3 


A033 

-0.1342 

0.1977 

0.1497 

-0.9707 

0 

0 

-4.6784 

-2.1477 

-0.0038 

0.8306 

0 

0 

0 

1.0000 

0 

0.1950 

0 

0 

0.5226 

-0.1951 

-0.0012 

-0.1905 

0 

0 

0 

0 

0 

1.0188 

-0.0001 

0 

3.6819 

0 

-0.7044 

0 

3.6819 

-0.0002 

B033 

0.0005 

0.0432 

0.1342 

0.9542 

0.7603 

4.6784 

0 

0 

0 

0.0855 

-0.6293 

-0.5226 

0 

0 

0 

0 

0 

0 

C033 

-0.1342 

0.1977 

0.1497 

-0.9707 

0 

0 

-0.0154 

0.0007 

0 

0.0012 

0 

0 

0.5226 

-0.1951 

-0.0012 

-0.1905 

0 

0 

1.0000 

0 

-0.1913 

0 

1.0000 

0 

D033 

0.0005 

0.0432 

0.1342 

0.0001 

0.0049 

0.0154 

0.0855 

-0.6293 

-0.5226 

0 

0 

0 

A036 

-0.2137 

0.0309 

0.0907 

-0.9891 

0 

0 

-9.4091 

-3.7688 

-0.0108 

0.8289 

0 

0 

0 

1.0000 

0 

0.0311 

0 

0 

2.1602 

-0.1993 

-0.0057 

-0.2476 

0 

0 

0 

0 

0 

1.0005 

-0.0001 

0 

6.1857 

0 

-0.1923 

0 

6.1857 

-0.0002 

B036 

0.0019 

0.0686 

0.2137 

2.3575 

2.3503 

9.4091 

0 

0 

0 

0.1523 

-1.6656 

-2.1602 

0 

0 

0 

0 

0 

0 

C036 

-0.2137 

0.0309 

0.0907 

-0.9891 

0 

0 

-0.0411 

-0.0000 

0 

0.0020 

0 

0 

2.1602 

-0.1993 

-0.0057 

-0.2476 

0 

0 

1.0000 

0 

-0.0311 

0 

1.0000 

0 

D036 

0.0019 

0.0686 

0.2137 

0.0004 

0.0132 

0.0411 

0.1523 

-1.6656 

-2.1602 

0 

0 

0 

r 


I 

I 


A049 

-0.1355 

0.1866 

0.1537 

-0.9737 

0 

0 

-4.5050 

-2.3097 

-0.0039 

0.7497 

0 

0 

0 

1.0000 

0 

0.1826 

0 

0 

0.4975 

-0.2190 

-0.0008 

-0.1918 

0 

0 

0 

0 

0 

1.0165 

-0.0001 

0 

3.5935 

0 

-0.6455 

0 

3.5935 

-0.0002 

B049 

0.0003 

0.0422 

0.1355 

0.9720 

0.7344 

4.5050 

0 

0 

0 

0.0896 

-0.6034 

-0.4975 

0  0  0 
0  0  0 


C049 

-0.1355 

0.1866 

0.1537 

-0.9737 

-0.0151 

0.0008 

0 

0.0011 

0.4975 

-0.2190 

-0.0008 

-0.1918 

1.0000 

0 

-0.1796 

0 

D049 

0.0003 

0.0422 

0.1355 

0.0000 

0.0047 

0.0151 

0.0896 

-0.6034 

-0.4975 

0 

0 

0 

A052 

-0.2139 

0.0232 

0.0977 

-0.9899 

-8.1563 

-3.7618 

-0.0110 

0.7173 

0 

1.0000 

0 

0.0225 

1.9229 

-0.2004 

-0.0051 

-0.2365 

0 

0 

0 

1.0003 

5.7438 

0 

-0.1292 

0 

B052 

0.0014 

0.0648 

0.2139 

1.9950 

2.0730 

8.1563 

0 

0 

0 

0.1338 

-1.4570 

-1.9229 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1.0000 


0 

0 

0 

0 

-0.0001 

5.7438 


0 

0 

0 

0 


-0.0002 


C052 

-0.2139 

0.0232 

0.0977 

-0.9899 

0 

0 

-0.0382 

0.0001 

0 

0.0018 

0 

0 

1.9229 

-0.2004 

-0.0051 

-0.2365 

0 

0 

1.0000 

0 

-0.0225 

0 

1.0000 

0 

D052 

0.0014 

0.0648 

0.2139 

0.0002 

0.0116 

0.0382 

0.1338 

-1.4570 

-1.9229 

0 

0 

0 

A-5 


o  o  o  o  o 


A065 

-0.1308 

0.1779 

0.1592 

-0.9768 

0 

-4.5147 

-2.1408 

-0.0023 

0.9537 

0 

0 

1.0000 

0 

0.1723 

0 

0.5158 

-0.2069 

-0.0011 

-0.1594 

0 

0 

0 

0 

1.0147 

-0.0001 

3.4757 

0 

-0.5902 

0 

3.4757 

B065 

0.0002 

0.0409 

0.1308 

0.9389 

0.6949 

4.5147 

0 

0 

0 

0.0867 

-0.5673 

-0.5158 

0 

0 

0 

0 

0 

0 

C065 

-0.1308 

0.1779 

0.1592 

-0.9768 

0 

-0.0141 

0.0009 

0 

0.0009 

0 

0.5158 

-0.2069 

-0.0011 

-0.1594 

0 

1.0000 

0 

-0.1698 

0 

1.0000 

D065 

0.0002 

0.0409 

0.1308 

0.0000 

0.0044 

0.0141 

0.0867 

-0.5673 

-0.5158 

0 

0 

0 

A068 

-0.2093 

0.0008 

0.1003 

-0.9910 

0 

-7.4918 

-3.4403 

-0.0052 

0.9763 

0 

0 

1.0000 

0 

-0.0013 

0 

2.0306 

-0.1696 

-0.0051 

-0.2089 

0 

0 

0 

0 

1.0000 

-0.0001 

5.5965 

0 

0.0072 

0 

5.5965 

B068 

0.0013 

0.0633 

0.2093 

2.0111 

2.0120 

7.4918 

0 

0 

0 

0.1304 

-1.3927 

-2.0306 

0 

0 

0 

0 

0 

0 

-0.000 


-0.0002 


C068 

-0.2093 

0.0008 

0.1003 

-0.9910 

0 

0 

-0.0364 

0.0004 

0 

0.0016 

0 

0 

2.0306 

-0.1696 

-0.0051 

-0.2089 

0 

0 

1.0000 

0 

0.0013 

0 

1.0000 

0 

D068 

0.0013 

0.0633 

0.2093 

0.0002 

0.0110 

0.0364 

0.1304 

-1.3927 

-2.0306 

0 

0 

0 

A-6 


oooootsi  oooo  ooooo 


A081 

-0.1441 

0.1484 

0.1627 

-0.9830 

0 

0 

-4.6721 

-2.1018 

-0.0008 

0.9807 

0 

0 

0 

1.0000 

0 

0.1403 

0 

0 

0.4422 

-0.2068 

-0.0023 

-0.1538 

0 

0 

0 

0 

0 

1.0098 

-0.0001 

0 

3.4168 

0 

-0.4748 

0 

3.4168 

-0.0002 

B081 

0.0002 

0.0402 

0.1441 

0.9445 

0.6894 

4.6721 

0 

0 

0 

0.0876 

-0.5491 

-0.4422 

0 

0 

0 

0 

0 

0 

C081 

-0.1441 

0.1484 

0.1627 

-0.9830 

0 

0 

-0.0153 

0.0010 

0 

0.0008 

0 

0 

0.4422 

-0.2068 

-0.0023 

-0.1538 

0 

0 

1.0000 

0 

-0.1390 

0 

1.0000 

0 

D081 

0.0002 

0.0402 

0.1441 

0.0000 

0.0043 

0.0153 

0.0876 

-0.5491 

-0.4422 

0 

0 

0 

A084 

-0.2350 

-0.0337 

0.1030 

-0.9912 

0 

0 

-8.0207 

-3.2912 

-0.0054 

0.9187 

0 

0 

0 

1.0000 

0 

-0.0352 

0 

0 

1.8108 

-0.1275 

-0.0045 

-0.2228 

0 

0 

0 

0 

0 

1.0006 

-0.0001 

0 

5.4492 

0 

0.1917 

0 

5.4492 

-0.0002 

B084 

0.0014 

0.0618 

0.2350 

1.5471 

1.9663 

8.0207 

0 

0 

0 

0.1094 

-1.3270 

-1.8108 

0 

0 

0 

0 

0 

0 

C084 

-0.2350 

-0.0337 

0.1030 

-0.9912 

0 

0 

-0.0398 

0.0002 

0 

0.0014 

0 

0 

1.8108 

-0.1275 

-0.0045 

-0.2228 

0 

0 

1.0000 

0 

0.0352 

0 

1.0000 

0 

D084 

0.0014 

0.0618 

0.2350 

0.0002 

0.0105 

0.0398 

0.1094 

-1.3270 

-1.8108 

0 

0 

0 

A-7 


A097 

-0.1413 

0.1216 

-4.4734 

-2.0836 

0 

1.0000 

0.5101 

-0.2114 

0 

0 

3.3284 

0 

B097 

-0.0001 

0.0392 

0.9606 

0.6702 

0 

0 

0.0904 

-0.5217 

0 

0 

0 

0 

C097 

-0.1413 

0.1216 

-0.0146 

0.0012 

0.5101 

-0.2114 

1.0000 

0 

D097 

-0.0001 

0.0392 

-0.0000 

0.0041 

0.0904 

-0.5217 

0 

0 

AlOO 

-0.2165 

-0.0394 

-7.6192 

-2.9495 

0 

1.0000 

1.5893 

-0.1280 

0 

0 

5.0075 

0 

BlOO 

0.0009 

0.0573 

2.0303 

1.6786 

0 

0 

0.1310 

-1.1371 

0 

0 

0 

0 

ClOO 

-0.2165 

-0.0394 

-0.0337 

0.0006 

1.5893 

-0.1280 

1.0000 

0 

DlOO 

0.0009 

0.0573 

0.0001 

0.0089 

0.1310 

-1.1371 

0 

0 

0.1677 

-0.9892 

0 

-0.0005 

1.0067 

0 

0 

0.1109 

0 

-0.0042 

-0.1307 

0 

0 

1.0061 

-0.0001 

-0.3669 

0 

3.3284 

0.1413 

4.4734 

0 

-0.5101 

0 

0 


0.1677 

-0.9892 

0 

0 

0.0005 

0 

-0.0042 

-0.1307 

0 

-0.1102 

0 

1.0000 

0.1413 

0.0146 

-0.5101 

0 

0.1120 

-0.9926 

0 

-0.0052 

0.9048 

0 

0 

-0.0434 

0 

-0.0053 

-0.1806 

0 

0 

1.0009 

-0.0001 

0.2169 

0 

5.0075 

0.2165 

7.6192 

0 

-1.5893 

0 

0 


0.1120 

-0.9926 

0 

0 

0.0010 

0 

-0.0053 

-0.1806 

0 

0.0433 

0 

1.0000 

0.2165 

0.0337 

-1.5893 

0 


-0.000 


-0.000 


A-8 


OOOO  KJOOOOO  oooo  roooooo 


Appendix  B:  MATLAB  User  Written  Routines 


%  routine  dovl 
%  builds  al,bl,cl,dl 

\  inputs:  wla4-wla5,wlb4-wlb5,wlk4-wlk5 
\  outputs:  al,bl,cl,dl 

\ 

%  W1  on  each  channels  4,5  Is  of  the  form  k*(S'i-b)/(s-t-a} 

\  W1  on  1-2  is  unity 

num4=vlk4*(l  vla4]; 

num5=vlk5*(l  vla5]; 

den4=[l  vlb4]; 

den5=[l  wlb5]; 

al=zeros (2,2); 

bl=zeros (2,4); 

cl=zeros (4,2); 

dl=zeros (4,4); 

dl(l,l)=l; 

dl(2,2)=l; 

Ia,b,c,d ]=tf 2ss (num4,den4 ) ; 

al(l,l)=a; 

bl(l,3)=b; 

cl(3/l)=c; 

dl(3,3)=d; 

clear  abed; 

(a,b,c,d]=tf 2ss (num5,den5) ; 

al(2,2)=a; 

bl(2,4)=b; 

cl(4,2)=c; 

dl(4,4)=d; 

clear  abed  nuro4  numS  den4  den5 


\  forms  weight  We  on  outputs  betadot , Ny, pslgt,  and  y. 
\  inputs:  wekl-wek5  veal-vea5  webl-web5  vecl&vec3 
%  outputs:  matrices  ae  be  ce 
\  routine  dove 

%  We  for  betadot  and  ny  is  We  =  k*a/(s+b) 

%  We  for  psigt  and  y  is  We  =  k*(s+a)/(s+b) (s+c) 

\  weight  on  betadot 
num=vekl*veal; 
dens[l  vebll; 

(a,b,c,d ]«tf 2ss (num^den) ; 

ae (1, 1 )=a; 

be(l,l)=b; 

ce(l,l)=c; 

clear  num  den  a  b  c 

%  weight  on  Ny 

num=vek2  *  vea2; 

den=[l  veb21; 

[a,b^c,d]-tf 2ss (num^den) ; 

ae(2,2)=a; 

be(2,2)=b; 

ce(2,2)=c; 

clear  num  den  a  b  c 

%  weight  on  psigt 

num=  vek4*(l  wea41; 

den=[l  (web4+wec4)  web4*wec4]; 

(a,b,c,d]=tf2ss(num,den) ; 

ae(3:4,3:4)=a; 

be(3:4,3)=b; 

ce ( 3, 3 : 4 )=c; 

clear  num  den  a  b  c 

%  weight  on  y 

num=  wek5*ll  weaSl; 

den=[l  (web5+wec5)  web5*wec5I; 

[a,b,c,d  ]=tf 2ss (num,den) ; 

ae(5:6,5:6)=a; 

be(5:6,4)=b; 

ce(4,5:6)=c; 

clear  num  den  a  b  c 


B-2 


%  routine  dove 

%  forms  weight  We  on  outputs  betadot,Ny,psigt,  and  y. 
\  inputs:  wekl-wek5  veal-veaS  webl-web5  vecltvec3 
%  outputs:  matrices  ae  be  ce 
\  We  for  betadot  and  ny  is  We  »  k*a/(s+b) 

\  We  for  psigt  and  y  is  We  *  k*(s+a)/(s+b) (s+c) 

\  weight  on  betadot 
num«vekl*veal; 
den-[l  vebl]; 

Ia,b,c,d]«t£2ss(num/den) ; 

ae(l,l)=a; 

be(l,l)=b; 

ce(l,l)*c; 

clear  num  den  a  b  c 

\  weight  on  Ny 

num=vek2  *  vea2; 

den-[l  veb2]; 

[a,b,c,d ]=t£2ss (num,den ) ; 

ae(2,2)=a; 

be(2,2)=b; 

ce(2,2)*c; 

clear  num  den  a  b  c 

\  weight  on  psigt 

num=  vek4*[l  vea4]; 

dens(l  (web4+wec4)  web4*wec41; 

(a,b,c,d]st£2ss(num,den) ; 

ae(3:4/3:4)~a/ 

be(3:4,3)=b; 

ce ( 3/ 3 : 4 ) =c; 

clear  num  den  a  b  c 

%  weight  on  y 

num=  vek5*[l  vea5]; 

den-ll  (web5+vec5)  web5*wec5]; 

(a,b,c,d ]=t£2ss (num,den ) ; 

ae(5:6,5:6)=a; 

be (5:6,4 )=b; 

ce(4,5;6)=c; 

clear  num  den  a  b  c 


%  routine  doabcd 

\  builds  k  B1  B2  Cl  C2  Dll  D12  D21  D22  natrices 

\  nsa  -  i  states  in  the  plant 

%  nsl  -  I  states  in  VI 

\  nse  -  I  states  in  VE 

\  ndl,nd2,nd3  -  I  exogenous  inputs 

%  nu  -  I  o£  control  inputs  to  plant 

\  ne  -  I  of  controlled  inputs 

\  nog  -  I  of  plant  outputs 

%  rhoda^rhodr  -  penalties  on  control  inputs 

nsa-8; 

nsl»2; 

nse=6; 

ndls4; 

nd2=l; 

nd3®4; 

nu=2; 

nel=4; 

ne2=2; 

nog=4; 

rhodadrs[rhoda  0;0  rhodrl; 

\  Bake  A 

A  =  [a084  zeros(nsa,nsl)  zeros (nsa, nse ) ; 

zeros(nsl,nsa )  al  zeros (nsl, nse } ; 

-be*c084  be*cl  ae  J; 

%  make  B1 

B1  s  [zeros(nsa,ndl)  bg084  zeros (nsa, nd3 ) ; 

bl  zeros(nsl,nd2)  zeros (nsl, nd3 ) ; 

be*dl  -be*dg084  -be*d3  ]; 

\  make  B2 

B2  »  [b084  ; 

zeros(nsl  nu)  ; 
zeros(nse,nu} ] ; 

%  Bake  Cl 

Cl  =  ( zeros (nel, nsa )  zeros (nel, nsl )  ce  ; 

zeros(ne2,nsa)  zeros (ne2, nsl )  zeros ( ne2, nse )] ; 
\  Bake  Dll 

Dll  s  [zeros(nel,ndl)  zeros (nel, nd2 )  zeros(nel,nd3); 

zeros (ne2,ndl )  zeros(ne2,nd2)  zeros (ne2,nd3) ) ; 

\  Bake  D12 

D12  s  Izeros(nel,nu)  ; 

rhodadr*eye(2) 1; 

%  Bake  C2 

C2  «  [-C084  cl  zeros (nog, nse) ]; 

%  Bake  D21 

D21  «  [dl  -dg084  -d31; 

%  Bake  D22 

D22  s  Izeros(nog,nu) ] ; 
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%  routine  close 

\  closes  loop  with  plant  and  compensator 

\  ndl,nd2,nd3  -  Inputs  to  closed  loop  system 

\  no  -  outputs  of  closed  loop  system 

\  nsc  ”  i  of  states  in  the  compensator 

\  a,b,c,d  -  closed  loop  system  matrices 

ndl«4; 

nd2=l; 

nd3*4; 

no*4; 

nsc=16; 

%  make  a 

a  =  [  a084  b084*ccp; 

-bcp*c084  acp  1; 

\  make  b 

b  «  {zeros(nsa«ndl}  bg084  zeros (nsa,nd3) 

bcp  -bcp*dg084  -bcp*d3  ] 

%  make  c 

c  -  [c084  zeros (nO/Hsc) ] ; 

\  make  d 

d  =  [zeros(no,ndl)  dg084  d31; 


%  routine  zh2064 

%  builds  matrices  for  computation  of  H2  compensator 

\  calculates  compensator,  close  loop,  and  runs  simulation 

\  plots  states  and  outputs  of  system 

dove; 

dovl; 

abcd064; 

h21qg; 

close084; 

cleanup 

pack 

[y,x]slsim(a,b,c,d,u,t,xO) ;x=x( :,1:8) ; 

plot(t,x( :,1) );grid;ylabel( 'beta* );xlabel( 'time  -  sec');pause 
plot(t,x( : ,2) );grid;ylabel( 'p' );xlabel( 'time  -  sec');pause 
plot(t,x( :,3) );grid;ylabel( 'phi* };xlabel( 'time  -  sec');pause 
plot(t,x( :,4) );grid;ylabel( 'r ' );xlabel( 'time  -  sec');pause 
plot (t,x( : , 5} } ;gr id;ylabel ( *psi ' ) ;xlabel { ' time  -  sec');pause 
plot(t,x( :,6) );grid;ylabel( 'y' );xlabel( 'time  -  sec');pause 
plot(t,x( :,7) );grid;ylabel( *dela* );xlabel( 'time  -  8ec'};pause 
plot(t,x( :,8) );grid;ylabel( 'delr ' );xlabel( 'time  >  sec');pause 
plot ( t,y( :, 1 )) ;grid;ylabel ( 'beta  dot ') ;xlabel (' time  -  sec'};pause 
plot(t,y( :,2) };grid;ylabel( 'ny' );xlabel( 'time  -  sec');pause 
plot(t,y( :,3) ); grid; y label ( 'psigt ' ) ;xlabel ( ' time  -  sec* ); pause 
plot ( t,y( : , 4 ) ) ;gr id;ylabel ( 'y ' ) ;xlabel ( ' time  -  sec');pau8e 
save  xytemp  x  y  a  b  c  d  acp  bcp  ccp  dcp  t 
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%  routine  pltvlp 

%  plots  the  bode  plot  of  the  input  weight  on  psigt 
nua*vlk4*ll  vla4];den«[l  vlb41;pltbode;gzid; 
ylabelCnag  VI (PSIGT)  -  db') 
xlabeK’FREQ  -  rad/sec') 


\  routine  pltbode 

\  plots  the  Bode  plot  of  a  transfer  function 

(Bag,p]sbode(nu]&,den,  freq)  ;Bag-20*logl0  (mag) ;  semi  logx(  freq^aag) 
grid 


\  routine  svplts 

%  plots  the  loop  gain  singular  value  plot,  the  sensitivity 
%  function  singular  value  plot,  and  the  complementary 
%  sensitivity  singular  value  plot 

[a,b,c,d]«ser ies (acp,bcp,ccp,dcp,a084,b064,c084,d084 ) ; 

svg]cssigma(a,b,c,d,l,v) ; 

svg)cs20*logl0(svg)( ) ; 

svg)(«svg)c(l:2, ; ); 

svs=sigma (a,b,c,d, 3,v) ; 

svs=-20*logio (svs ) ; 

(a,b,c,d ]«feedbk (a,b,c,d, 2 ) ; 

svcl-sigma(a,b,c,d,l,v)  ,* 

svcl=20*logl0 (svcl ) ; 

svcl=svcl(l:2, : ); 

semi logx ( V, svgk ); grid; 

titleCLoop  Gain  SV');pause; 

semi logx (V, svs ); grid; 

title ( 'Sensitivity  Func  SV');pause; 

semi logx (v, SVC 1 ) ;gr id; 

title ( 'Closed  Loop  SV) 

save  svtemp  v  svgk  svcl  svs 
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